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MSTRACT 


In Rosen's "blmetric" theory of gravity the (local) speed of gravi- 
tational radiation v^ is determined by the combined effects of cosmolog- 
ical boundary values and nearby concentrations of matter. It is possible 

for v to be less than the speed of light. I show here that emission of 
S 

gravitational radiation prevents particles of nonzero rest mass from 
exceeding the speed of gravitational radiation. Observations of relativ- 
istic particles place limits on v^ and the cosmological boundary values 
todav, and observations of synchrotron radiation from compact radio 
sources place limits on the cosmological boundary values in the past. 
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LIST OK SYMBOLS 
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V 

H 
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e 

P 

6 
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Greek lower case mu 
nu 

" " •' eta 

alpha 

ganina 

" curly epsilon 

beta 

delta 

•' " " lambda 

" " " pi 


V = nabla 

Q * Greek lower case theta 

= integral 

5 = upper case script S 

u: = Greek lower case omega 

( a lower case script "ell" 
A = Greek upper case delta 

OO S3 infinity 

W = Greek upper case omega 

? a " lower case xi 
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1. INTRODUCTION 


Several years ago Nathan Rosen (1] proposed a now theory of gravity, 
the "blmotrlc" theory — the two metrics being the physical metric and 
a flat, "background" metric n . 1 The theory Is perhaps better described 
as a two-tensor, metric theory (see [2| for discussion). It Is a metric 
theory in the sense that the physical metric obeys the Einstein Equivalen e 
Principle: in the local, freely tailing frames of g the nongravitut lonel 

laws of physics reduce to those of special relativity. One immediate con- 
sequence is local conservation of nongravi tat ional stress-energy T — the 
matter-response oquat ion: 


T llV - 0 (11) 

;v 

The auxiliary, symmetric cwo-tensor n v can he thought of as a second metric; it 
is constrained to be flat and is used in constructing the field equations for 
the physical metric (see Section 215). In a series ot papers Rosen and others 
have analyzed various consequences of the theory, Including the maximum mass 
of neutron stars |3], cosmological models l 4 | , equations of motion (51, 
gravitational radiation |b], and other topics |7|. 

The traditional testing ground for such a theory Is the solar systei, 
where observations at today's accuracies probe the theory's predictions t.' 
pos t-Newton lan order (see (8) for a review). Lee et a_l. (^] have calculated 
the post-Newtonian limit of P.os»*n's theory and shown that it Is the same as 
that of general relativity, except for the preferred- frame PPN parameter 
| 'i a discussion of the Nordtvedt-Wi 11 Parametrized Post-Newtonian (PPN) 
Formalism and a description of the meaning of the PPN parameters, see charter 
w or | io | ; i n part 1 cu Lar, ox 39.5, >•. < * ton! 
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gravitational "conatant" G arc determined by the distant matter in the 
Universe, which reaches Into the solar system through boundary conditions 
applied far outside it. An appropriate adjustment of the cosmological 
boundary values brings the theory into agreement with present limits on 
iV, and on the time rate of change of G. Put the other wuy around, these 
limits place constraints on the possible boundary values. One wuy to tes' 
the viability of the theory is to construct cosmological models and ask 
whether the models can be made consistent with these constraints. In till • 
paper 1 point out a new set of observations which yield particularly 
stringent constraints on the cosmological models in Rosen's theory. 

The two metrics in Rosen's theory play different roles, Gravltatlor.il 
radiation propagates along "light" cones of the flat metric, while light 
propagates along "light" cones ot the physical metric. The two "light" 
cones need not coincide, so the speed of gravitational radiation is, in 
general, different from the speed of light. Lee et al. [0] showed that tae 
speed of gravitational radiation, as measured by an observer at rest in the 
universal rest frame far irom anv local concentration of matter, is dctoi - 
mined solely by the cosmological boundary values. Tills speed v ^ is 
related to >. 0 by 


C. 

V 

RC 


(l 



( 1.2a) 


!n the vicinity of a local source of gravity with (dimensionless) Newtonian 
potential U « l (U > 0) , the speed of gravitational radiation increases 
to 


v 

R 


v (l r 2U) 
8 C 


(1.2b) 


(sec Section :’W). It is possible for v to be less than the speed of llgnt. 


I show here (Section 311) that aw a particle of nonzero rest mass Is 
accelerated through the gravitational "light" cone. It emits an infinite 
amount of energy In gravlt atlonal radiation. It lollows that, if < 1, 
the speed of gravitational radiation is the ultimate speed for such par- 
ticles; they cannot escape the gravitational "light" cone. As a result, 
observation of a relativistic particle with Lorentz factor Y provides a 


lower bound for v at the point of observation: 
R 

•-v-K 2 • 


(1.1a) 


If the Newtonian potential at the point of observation Is known, one also 

obtains a lower bound for v . This lower bound cun be re-expressed as an 

gc 

upper bound on the value of ot ~ 2(1 - v ): 
rr 2 gc 


« 2 < Y ‘ + 4U 


(1.1b) 

Equations (1.3) are the basis for obtaining observational constraints on v 


and a, (see Section 4). 

In this paper l analyze the gravitational radiation emitted by par- 
ticles moving at speeds near the speed of gravitational radiation. This 
analysis leads to the conclusion that. In Rosen's theory, particles of 
nonzero rest mass cannot exceed the speed of gravitational radiation. This 
conclusion is likely to have far wider applicability than just to Rosen's 
theory. The detailed analysis presented here does not depend critically on 
any special feature of Rosen's theory; one can make a strong case that a 
similar anal) Is holds in any theory of gravity which permits the speed o f 
gravitational radiation to differ from the speed of light (see (11] for a 
brief review of such theories). Indeed, It seeira likely that the "gravita- 
tional speed limit" la a feature of ill such theories. 


Another crucial test of Rosen's theory comes from observations of the 


change In orbital period of the binary pulsar [12]. Unless the two compon- 
ents of the binary system have Identical ratios of gravitational binding 
energy to inertial mass, Rosen's theory predicts that the svstem will emit 
dipole gravitational radiation and th.it the radiation will carry away 
negative energy [13]. Observations of the binary pulsar are now good 
enough that Rosen's theory can be ruled out — unless the two ratios are 
the same to within less than a percent [14], 

Section 2 develops the formalism for analyzing gravl tat lonal radia- 
tion emission from weak-field systems In Rosen's theory. Section » hull, s 
upon this foundation to Justify the claim that a particle of finite rest 
mass cannot ex' _ed the speed of gravitational radiation. Section 3A calcu- 
lates the energy spectrum of gravitational Cherenkov radiation emitted by 

a particle moving with uniform veil city v > v , and Section 3B analyzes the 

R 

energy emitted as a particle is accelerated through the gravitational "light" 

cone. The result of these considerations is Kqs. (1.3), which Section 4 

uses to obtain observational limits on v and a~ (v ). Section 5 argues 

8 - gc 

that these constraints apply to any theory of gravity with a variable speed 
of gravitational rudlation. 


2. FOUNDATION FOR ANALYZING EMISSION OF GRAVITATI0NA1 i >1 ATT ON 

This section lays the foundation for analyzing emission of gravita- 
tional radiation from weak-field, linearized svstems in Rosen’s theory. The 
foundatior will be laid In two pieces: the first piece is construction of 

coordinates which take into account matching to boundary values provided Sv 
an external grav i t at ional field; the second piece is construct ion of equa- 
tions governing generation of gravitational radiation and specifying the 
amount of energy the radiation carries. 
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A . L s o 1 jtetl Sources and Prefer red Coordinates 






Below I shall ileal with "isolated" sources of gravity, such as the 
solar system or one of the ult rare lat 1 vis t ic particles of Section 3. Since 
such sources are not actually alone in the Universe, it is necessary to 
describe briefly what is meant by an "Isolated" source. 

The key feature of an isolated source Is that the gravitational field 
can be split into two pieces: the field of the Isolated source (the local 

field), which applies near the source; and the field of the reRt of the ir »t— 
ter in the Universe (the external field), which applies far from the source. 

To understand the conditions necessary for such a split, consider the 
length scales characteristic of the source and the external field. The si irce 
is characterized by two lengths: Its physical size R and the length Cm cjr- 
responding to its mass m. The external field is also character ized by tv > 
lengths: a typical radius of curvature a and the length L over which the 

external field varies appreciably. Let r be the distance from the sourr • at 
which the curvature produced by the source becomes comparable to the extt nal 
curvature: 

r = (Cn.a 2 ) 1/3 . (2.1) 

o 

To get a clean split between the local and external fields, the source must 

be buried deep inside r (R «.*. r ), and r must be much smaller than the »x- 
r o o o 

temal scales (r << mln{a,L}). These two conditions translate Into 
o 

Gm/R 3 » a' 2 , (2 2a) 

Cm << mini a, L(L/a)" 1 . (2.2b) 

The curvature produced by the isolated source is a large, but small-scale 
"bump" in the large-scale external curvature. 
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When Kqs. (2.2) are satisfied, the region around the source can be 
broken up Into three parts, which provide a natural split of the gravita- 
tional field: 

(I) the local-field region , In which the curvature of the Isolated 
source dominates: 

r < r , (2.3a) 

- o 

where r is distance from the source; 

(II) the transition region , In which the external curvature dominates, 
but which is small enough that the external field is nearly homogeneous: 

r < r < r. t • tnln(a,L) , (2. lb) 

o l 

where e is a suitably chosen factor less than one; 

(III) the external-field region : 

r > r^ . (2 , : c) 

The nearly fiat transition region splits the gravitational field Into local 

and external fields. The only connection between the two fields is the re- 
quirement that they match smoothly In the transition region; from the point 
of view of the local field, the external field establishes boundary condi- 
tions In the transition regLon. 

The boundary conditions are made explicit by choosinR a specific 
coordinate system. A particularly convenient set of coordinates can be 
constructed as follows. Consider the external gravitational field in the 
absence of the Isolated source. Let an observer falling freely in this 
field construct Fermi normal coordinates lx'} t |t,xM In the vicinity of 
his world line [13]. In these coordinates the two metrics can be expanded 
about the observer's world line: 


«*3 " *ad 


4 [terns of order K 


x^x k 'v (r/a)") 4 ••• , 


(2.4a) 


ajtfk 

n a ^(t,x^) - n^(t.x^-O) 4 n a|; k (t > x J .<»x k 4 ••• , (2.4b) 

where lx Miukowsklan, l.t M I! !' " dlag(-l # 4l,4l ,4 1) , and R ( J> . 

is the Rienuinn tensor derived from the physical metric. Now introduce the 
isolated source in the vicinity of the fiducial world line and use these 
coordinates to solve for its local field. The flat metric retains the form 
(2.4b) » and tha phyaical aatrl< retains the form (2.4a) in the transition 
region outside the source. 

Equations (2.4) display explicitly the boundary conditions to be applied 
in the transition region. In general relativity, which has only a physical 
metric, the external field influences the Isolated source only through the 
Riemaun and higher-order terms in Eq. (2.4a), which represent tidal and 
higher-order multipole forces on the isolated source. The situation Is dif- 
ferent in Rosen's theory because of the presence of the fl it metric. Although 
the region around the source has been split cleanly into local and extern:! 
parts, n^, cannot be adjusted independently In the two regions; rather, the 
external field determines the form of , in the transition region, and a 
particular choice of coordinates together with its flatness then determines 
ti In the local-field regie: (The above choice of coordinates Insures that 

n is nearly constant in the local field region.) In general, the external 
field prohibits finding coordinates such that both g and p ^ are nearly 
Miukowsklan In the transition region. This lack of "meshing" allows the ex- 
ternal field to reach into the vicinity of the Isolated source and affect 
local gravitation physics. (Will [2, Section 5.3! gives a general discussion 
of tiie manner in which auxiliary tensor fields In metric theories of gravity 
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0 


couple local gravitation physics to an external field. ) 

Gravitational radiation emitted by the source la analyzed In the tran- 
sition region. In order to separate the radiation from the external curva- 
ture, the wavelength \ of the radiation must be much smaller than external 
scales : 

X << mln(a,L} , (2.5) 

a requirement which also guarantees that the wave zone of the radiation ex- 
tends into the transition region. In Section 3 I will be Interested In 
calculating gravitational radiation emission In the linear approximation. In 
this limit another consequence of (2.5) Is that. In Eqs . (2. A), one can 
ignore both the tidal terms in and the spatial and temporal derivatives 
of '1^; these terms cannot affect radiation at wavelengths much smaller 
than their own characteristic lengths. 

As a result, in calculating the gravitational radiation emitted by an 
isolated source in the linear approximation, one can always use coordinates 
with the following two properties: 

Property 1 . The physical metric g ( is asymptotically Minkowskian in 
the transition region far from the source. 

Property 2 . The flat metric , is a nearly constant matrix in the 
the locil-field and transition regions; its sI.jwIv changing values 
are de'ermined by the external field, and its temporal derivatives 
can be Ignored. 

1 shall refar to a coordinate system which satisfies these two properties ar 
a preferred coordinate system . Such coordinates are particularly useful for 
analyzing gravitational radiation emission: Property 1 insures that the 

coordinates provide a good reference frame for an observer in the transition 


a 


region monitoring the emitted radiation; and Property 2 Insures that the 
field equations and gravitational stress-energy assume a particularly simple 
form (see Section 2B). Properties 1 and 2 do not uniquely specify the coor- 
dinates; Instead, they specify a family of preferred coordinate systems, 
the members of which are related by arbitrary Lorentz transformations and 
translations. Throughout the following I shall use preferred coordinates. 

Now restrict attention to the sources considered In Sections 3 and 4 — 
ultrareletivistlc particles moving in a typical astrophvslcal environment. 

For such sources, the external field must include both the smoothed-out 
cosmological solution and the fields of nearby, large-scale density enhance- 
ments. A typical source might be a cosmic-ray proton near the Earth; then 
the nearby density enhancements include the Virgo cluster, the Local Group, 
the Galaxy, the solar system, and the Earth. To sufficient accuracy the 
gravitational fields of the large-scale density enhancements can be treated in 
the weak- fie Id, s low-motion approximation. If the Universe is homogeneous and 
isotropic (assumed henceforth), the solution for the full external field in- 
cluding the cosmological boundary values and nearby density enhancements - is 
that given in reference [9], In the universal rest »rame — the frame In which 
the cosmological fluid is at rest the two metrics are given bv 


*00 " 

-1 + 2U , 

C.ha) 

8 0J " 

0 

(2. hb) 

g Jk- 

6 Jk (l+2U) 

(2.6c) 

ihjl 

, -1 -1 -1 -1, 

- diag(-c 0 , c 1 .Cj ,c 1 ) , 

(2.6d) 

and Cj 

are determined by the cosmological solution. 

and U is the 


Newtonian potential due to those nearby density enhancements which produce a 
significant deviation from the o sinologi cal solution. The r^ component have 
been neglected, since they are much smaller than U for s^ow-motlon ources. 


The external field (2.6) la used to construct preferred coordinates 
appropriate for analyzing an isolated source. For manv purposes the roost 
convenient set of preferred coordinates Is obtained by using a (freely falling) 
fiducial observer who is Initially at rest in the universal rest frame. In the 
resulting preferred coordinates i* asymptotically Minkowski an In tht trar.: i 

tion region (Property l), and ; ( . is given by 

n Q0 * * c q 1(1 + 2U) . (2. 7a) 

n 0j “ 0 , (2.7b) 

“ tj 1 6 Jk (l - 2U) . (2.7c) 

where U (= constant) is evaluated in the vicinity of the isolated source. 

These preferred coordinates will be called the local universal rest (LUkF)j 
they will be used for all the calculations in Section 3. 


B . Llnear l ? . u Field Equations and Gravitational Stress-Energy 

It is not necessary to give the full, nonlinear Rosen field equations 
here; only the linearized version will be needed. For the full equation-, 
the reader is referred to the original papers of Rosen [1] and to [9]. 

For a weak-fleld oource, the physical metric in a preferred coordinate 
system Is nearly Minkowskian In both the local-field and transition regions. 

In the usual way, define the metric perturbation h to be the deviation of 

UV 

g from Minkowski an: 

uv 


UV 


(B) 4 . h 
g + h , 

uv uv 


8^ = diag (-1 ,+l ,+l ,+l) ; 


( 2 . 8 ) 


and let h be the trace- reversed metric perturbation: 

uv 


where h 


T- _ . 1 0 (B) . 

h = h - r g h , 

UV pv 2 °uv 

(B)uv (i The j n< jf ce8 0 f h 

UV UV 


( 2 . )) 


and h 


UV 


are raised and lowered 


10 


using g 


(B) 

Uv 


To linear order In the metric perturbation, the field equations In 
any preferred coordinate .system are given by 


n aS h - -16* (-g/-n) 1/2 c T 

Mv.of} 6 o uv 


(2.10a) 




Here T) is the Inverse of ~ u ,.i S and 1] are the determinants of g^. and , 
respectively; and Is a coupling constant with dimensions of the Newtonian 
gravitational constant (see [9]). In LURF coordinates, Eq. (2.10a) becomes 


V 2 h - (1/v 2 ) h nn - -16* GT 
MV g UV.OO M^ 

1 / 2 , 


(2. 10b) 


Here G ■ (c^c^) G q is the gravitational "constant" at the epoch of 
interest (as measured, e.g., by a Cavendish experiment performed far awav 

from any local density enhancements), and v — the speed of gravitational 

2 

radiation in the LURF — is given by Eqs. (1.2). Equation (1.2a) uses the 

2 -1 

results of (9] to relate v to ol [v * c./c _ ■ (1 + a„) ). 

gc 2 gc 10 2 

The linearized matter- response equations are obtained from Eq. ('1.1): 


- 0 

,v 


( 2 . 11 ) 


Just as in genera’ relativity or any other metric theory, gravitational ef- 
fects disappear from the matter-response equations at linear order; the 
linear approximation is valid only so long as the motion of the source is 
governed by nongravitational forces. 

To analyze gravitational radiation emitted by a source, one must ha 
able to cal ul ite the energy and momentum carried by the radiation. Rosen 
[1) has demons t rated the existence of a stress-energy complex 0 which is 
conserved with respect to the flat metric: 


11 


(2. 12a) 


0,-0 

u |v 


<-*'-«> 1/2 < T u V * t u V> 


(2.1?b) 


Tlu* tensor t ^ In Interpreted as the gravltatlona stress-energy; it is a 

quadratic expression In first derivatives of g witn respect to n (g 114 , ). 

ViV tIV M v | a 

To lowest order In the metric perturbation, it is given in any preferred 
coordinate system by 

v l / < .l/2/.avB l.v otf) r 1 r r I 

■ CT- < - n/ '“ ) l 4 u" -2 { |, n I'h • 


(2.13) 


. r* (B)iiVj- 

where h = g 


UV 


Equations (2.12) can be Integrated to obtain conservation laws for t e 
total 4-momentum. In a preferred coordinate system, the 4-momentum P of the 
source is defined by 

r 


a 


it n 


(2 .14) 


P transforms like a 4-vector under Lorentc transformations among the pre- 

a 

ferred coordinates, and its indices are raised and lowered using g^ . Now 
surround the source with a closed 2-surface S which lies In the transition 
region, and let n be the unit outward normal (with respect to g ' to S . 
The conservation law (2.12) relates the change of 4-momentum inside $ to a 
■ix of 4-momentum thiough S 5 


dP_ 

dt 


- f (T aJ ♦ t uJ 


j . 

A 


)n , dA 


( 2 . 1 ^) 


Here the first index of t^ ' has been raised using g^. 

The important quantity for calculating energy loss due to gravitational 
radiation is t* the energy flux in the radiation. Its form (to lowest 


12 


order In the met ric perturbat Ion) Is particularly simple In LURF coordinates: 


t oj 





(2.16) 


Note that the field equations (2.10a) and the gravitational stress-energy 
(2.13) are not invariant unuer infinitesimal coordinate (gauge) transforma- 
tions. This lack of invariance reflects the fact that a gauge transformation 

destroys Property 2 of the preferred coordinates, i.e., does not remain 

(B) 

constant on local scales. In such coordinates, R v i y ^ 0 and terms con- 

/ n \ / d \ 

tainlrg | and g^ . | ^ c appear in the linearized i'.eld equations and in 
the gravitational stress-energy 


3. GRAVITATIONAL RADI/TION AND THE GRAVITATIONAL SPEF.D LIMIT 

In Rosen's theory the speed of gravitational radiation Is determined by 
the combined effects of the cosmological gravitational field and the gravita- 
tional fields of nearby, local concentrations of matter [F.qs. (1.2)]. 

Although the latter always tend to Increase v^, the cosmological field can 

force v to be less than the speed of light. It is this case — v <1 — 
g g 

that I consider in this section; in particular, I investigate the gravita- 
tional radiation emitted bv particles moving at speeds near v^. 

The motivation for doing so is provided by an analogy with electromag- 
net .m. A charged particle, moving through 3 material medium at a speed 
faster than the speed of light in the medium, emits electromagnetic Cherenkov 
[16] radiation. In any real medium dispersion restricts the radiation to a 
finite rang* of frequencies; however, for an idealized, dispersionless merium, 
the energy emitted diverges. Similarly, in Rosen's theory, a particle which 
exceeds the speed of gravitational radiation ought to emit g ravltat lonal 
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Cherenkov radiation. Moreover, the gravitational "medium" JLa dispersion- 
less (at least at high frequencies), so the electromagnetic analogy suggests 
that the energy emitted ought to diverge. If so, tills result would suggest 
that particles cannot exceed the speed of gravitational radiation. 

These ideas were first considered In a different context by Alchelburg, 
Ecker, and Sexl [17]. They considered a particle whose equation of motion 
apparently allows It to exceed tin* speed of light, but which is coupled to 
a field that propagates at the speed n f light. They argued that radiation 
reaction prohibits accelerating the particle to speeds greater than the 
speed of light. They showed, for example, that if such a particle is 
charged, the electromagnetic radiation it omits diverges as it is acceler- 
ated through the light cone. The situation considered here is similar, ..nd 
the analysis is patterned after their work. 1 shall first consider the 
gravitational Cherenkov problem and then analyze the power radiated in 
gravitational radiation as a particle is accelerated through the gr i.ta 
tional "light" cone. 


A. Gravl tational Che ronkov Had i at Lon 

Consider a particle with rest mass in moving with uniform velocity v 

relative to the LURF; let v v . In the case of Interest, v is very 

g K 

close to the speed of light. Adopt LURF coordinates and solve for the grav'- 
tational field in the linear approximation. The solution of the field equa- 
tions (? . 10b) proceeds exactly as in the analogous electromagnetic problem 
(.see, e.g. . [18], Section 14.9). The metric perturbation h forms a shock 
front along a cone which extends back from the instantaneous position of the 
particle (see Fig. 1); the angle 0 between the velocity v and the normal to 
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vanishes ; 


the cone is given bv cos 0,. ■ (v /v). Outside the cone h 

t R pv 


inside the cone. 




8 Gym 

o 

| x -vt|[l - (v/v )^sln“o]^ 


(3.1a) 


ll ■ — li 

*0j 00 


(3.1b) 


u j k 

h ■ h V v 

Jk 00 


(3 Ic) 


where u is the angle between the observation point x and the velocity v, and 
■) -1/2 

y = (1 - v“) is the particle's Lorentr factor. The field (3.1) represents 

gravitational Cherenkov radiation propagating in the direction normal to the 
Cherenkov cone, lly evaluating the energy flux using (2.16) and then following 
the procedure used for electromagnetic Cherenkov radiation, one obtains the 
energy d“E radiated into an angular frequency interval dn as the particle 
moves a distance d4: 


ti** t 2 'i 

~ 7 ~r r ■ Cm“u)( v - ' - 1) for v > v 
divdi o ' g 


0 . 2 ) 


This expression is similar to the Frank-Tamm [191 result for electromagnetic 
Cherenkov radiation. 

Equation (3.2) does not, of course, hold for all frequencies, and it 
is important to determine its region of validity. In using the formalism 
of Section 2, the above analysis neglects variations In the external gravi- 
tational field. However, since the particle is assumed to radiate for at 
infinite amount of time, these variations cannot be ignored; their effect is 
to modify Eq. (3.2) at low frequencies. To estimate the frequency at which 
such modification becomes important, consider a particle which radiates tor 


1 . 


only a finite time T *>* (r^/v) |h«* Kq. (2.3b)]. Then the particle's motion 
and the radiation it emits can he analyzed within the transition region, 
where the formalism developed in Section 2 is applicable. The emitted radia- 
tion is a pulse which lies just inside the Cherenkov cone (see Fig. 1). It 
is easy to show that, when the radiation is analyzed at a distance ^ r^ from 

the purticle's trajectorv, the pulse has a duration At ^ [(v/v )- l](T/8). 

R 

Thus the energy spectrum will be given by Eq. (3.2) for frequencies 

w > u' c = (1/At). The wavelength corresponding to the critical frequency & c 

is 


X 'v* e(v 1 - v *) • mint a, L} » jc (Y ? -Y “) • mints, L) , (3.3) 

c 8 — R 

2 - 1/2 

where y = (l - v ) . I have confirmed this result bv a detailed analysis 

8 R 

of the radiation emitted by a particle which moves faster than v for onlv a 

R 

fit. ice time. 

Variations in the external gravitational field can be regarded as produc- 
ing "dispersion" in the propagation of gravitational radiation. This disper- 
sion modifies the Cherenkov spectrum at frequencies below Note that, as 

v approaches v^, u> increases and dispersion affects more of the spectrum. 

Even when Y is not close to Y . the critical wavelength X is typically rather 

8 c 

small. For example, for a particle near the earth, the relevant external 

9 

scale is of order the radius of the earth: L % 10 cm. Using the smallest 

value of v allowed by the limits obtained in Section A [see Eq. (A.2b)J 
8 

-2 -11 2 
and choosing e * 10 , one finds A^ % (10 cm) [ 1 - (y^/y ) 1 • 

The importance of the preceding analysis lies not so much in estimating 


the size of > , but rather in demonstrating that, as long as v > v >t then is 
c g 

a finite critical frequency above which the gravitational "medium" is 
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dispersionless and Eq. (3.2) applies. In the purely classical analysis 
given above, the validity of F.q . (3.2) extends to arbitrarily high frequen- 
cies, and the spectrum diverges as u> ■* However, quantum mechanics often 
eliminates classical divergences, and one might expect a proper quantum- 
mechanical treatment to modify the classical spectrum at very high frequen- 
cies. In particular, conservation of energy might seem to require that the 

spectrum be cut off at a frequency u> - (>m /ft) corresponding to emission 

max o 

of a graviton whose energy is equal to the particle's energy. Applying tnis 
cutoff to Eq. (3.2) (and assuming (l ’ n dX >> w ) , one finds an energv loss rate 


dE! 

de 


Villi . - 

~~r 10' 1 " eV 

2v^fi 


-1 

cm 


for protons . 


3.4) 


Tills energv loss rate is so small that, if there is a cutoff at ui , the 

max 

effects of gravitational Cherenkov radiation are negligible even on galactic 
distance scales. 

However, the existence of the cutoff is by no means certain. The 
uncertainty arises because it is not clear that Rosen's theory, even in 'ts 
linearized version, can be quantized; the linearized field equations (2.10a) 
are not those of a canonical field theory. The difficulties that thereby 
arise are perhaps most apparent in an examination of plane gravitational 
waves in Rosen's theory: 

(i.) The Riemann tensor derived from an arbitrary plane wave has six 
independent polarizations — the most general polarization structure all wed 
in a metric theory. Even in the case v * l, where the theory is Lorentz- 
invariant, these six polarizations form a nonunitary representation of the 
inhomogeneous Lorentz group; they cannot be associated with massless quaita 
of definite, Lorentz-invariant helicity (see [20) for a general discussion 
of these issues). 
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(li) The time-averaged energy density In an arbitrary plane wave 
(v t 4 1), evaluated using (2.13), can be regarded as a quadratic form In 
the amplitudes of the ten independent potentials h^. (None of these ten 
potentials can be removed from the energy density by a gauge t » ans format 1 on ; 
see Section 2B.) When this quadratic form is diagonalized, one finds 
that four of the eigenvalues are negative. In other words, four of the ten 
degrees of freedom in the wave carry negative energy. The presence of 
negative-energy radiation has been noted previously in analyses of radiation 
emitted by binary systems in Rosen's theory 113]. In a theory with such 
negative-energy radiation, the stability of the vacuum Is uncertsin. 

Any attempt to quantize Rosen's theory must confront these two prob- 
lems. Even if they can be overcome, the presence of negative-energy radia- 
tion removes the raison d'etre for a cutoff at <a . The classical Cher-nkov 

■ max 

radiation (3.1) is made up of both positive- and negative-energy parts, the 
total energy emitted being a balance between the two. Quantum-mechanical ly , 
this Cherenkov emission might well be represented by multi-graviton proc 
esses in which both positive- and negative-energy gravitons are emitted. In 
such processes, conservation of energy imposes no restrictions on the fr. - 
quency of the emitted gravitons. 

Another potential quantum-mechanical cutoff is the Planck frequency 
- 1/2 

ui = (Oft) . If the Cherenkov spectrum (3.2) is cut off at , the 
p P 

energy loss rate becomes 


dE 

dli 


ym 

1 o 

2y 3 (tf/m o ) 


Y m p 

3 . c ,-14 . 

y (5x10 cm) 


for protons . 


(3.5) 


Just as for the cutoff at w , It is not clear that this cutoff should be 

max 

Imposed. However, even if it is, the loss rate (3.5) is large enough that 
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the limits obtained in Section 4 are not affected. 


The classical analysis of gravitational Cherenkov radiation hints 
at a serious problem in Rosen's theory. The divergence of the spectrum as 
u> •* « means that the energy emitted is Infinite (and positive). This result 
strongly suggests that particles cannot exceed the speed of gravitational radi- 
ation. It is not clear that a quant tun-mechanical treatment will eliminate 
the divergence, nor Indeed that such a treatment can be given. Even as a 
purely classical analysis, the above calculation has serious difficulties: 
it Is clearly Inconsistent and. Just as clearly, the linear approximation 
is not valid. However, there is little point in trying to patch up these 
difficulties. If particles cannot exceed the speed of gravitational raola- 
tion, a consistent calculation of gravitational Cherenkov radiation is not 
possible. More realistic and more relevant would be an examination of wnat 
happens as a particle is accelerated up to the speed of gravitational radia- 
tion. 

Before turning to this problem, it is interesting to ask about the 
Cherenkov radiation emitted by photons and other zero rest-mass particles. 

The best that can be done using the above calculation is to model a free 
photon as the limit v + 1, y ■* -, ym ( ♦ constant. Applying this limit to F.q. 
(3.2), one finds that free photons apparently do not produce any gravita- 
tional Cherenkov radiation. 

B . Acceleration through the Gravitational "Light" Cone 

Now consider a particle with rest mass m which has velocity v(t) 
in LURE coordinates. The particle is being accelerated h\ interactions 
with other matter and nongravi tat lonal fields. The objective is to 
evaluate, in the linear approximation, the energy emitted In gravitational 
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radiation as v approaches v^. In doing so, one must remember that the total 
stress-energy is conserved [Kq. (2.11)]. This means that one cannot, In 
general, neglect the radiation emitted by the matter and nongravltatlonal 
fields as they "recoil" from the Interaction. However, I shall argue that 
in the case of Interest here, this "recoil" radiation can be neglected. 

Imagine the following scenario for accelerating the particle — a 
scenario similar to those often envisioned for accelerating cosmic rays [21]. 
The particle Is accelerated by a series of "collisions" with local concen- 
trations of stress-energy. These "blobs" of stress-energy have masses much 

larger than m , and their velocities — both center-of-mass and internal — 
o 

relative to the LURF are small. In each collision, the momentum exchanged 
is small compared to the particle's momentum. The subsequent motion of 
the "blob" occurs on time scales much longer than the collision time; 
clearly, the radiation emitted by the "blob" does not diverge. Now consider 
the final stage of the acceleration process, when after many collisions 
the particle has attained a velocity so close to v^ that one more collision 
can push its velocity above v^. From the point of view of the partlculni 
"blob" involved, this collision is no different from the preceding ones. 
However, the radiation produced by the particle in this collision is beamed 
in the direction of its velocity, and the radiation diverges in that direc- 
tion as v approaches v . Therefore, in analyzing the final stage of the 

6 

acceleration process, one can neglect the "recoil" motion and calculate 

the radiation emitted by considering only the particle's motion. The 

results obtained will be valid when v is very close to v . 

8 

The field equations (2.10b) for a single-particle source can be 
solved in the same way as in electromagnetism (see, e.g., [18], Sec. 14.1). 
The h^ have the same form as the Lienard-Wiechert potentials. The energy 
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flux is evaluated using Eq. (2.16), and an Integral over a sphere In the 
transition region gives the power radiated: 


-4 nr 6 (v» v) 2 (-11 ♦ 26u ? - 9uS ♦ y 4 v 2 ( 1 - u 2 )(-11+ 12u 2 )J 
°-o v g[^ 6( yJ )2 < 1 - li2) " 4 + 3 Y 4 v 2 (i- U 2 )‘ 3 J 


where M = Y/Y , and where the last expression contains the leading-order 
8 

terms in the limit v ■* v . As anticipated, the power radiated diverges. 

8 

In a real situation, of course, the power radiated cannot diverge; Instead, 
radiation reaction diverges and prevents the particle from exceeding the 
speed of gravitational radiation. 


This calculation suffers from some of the same difficulties as th*- 


Cherenkov calculation. The particle radiates substantial amounts of energy 
only at very high frequencies where quantum corrections might well be Im- 
portant. For the reasons given earlier, the effect of these corrections 
is uncertain, and I shall ignore them. A perhaps more serious objection 
is that the linear approximation Is not valid; however, it seems unllkejy 
that the nonlinear terms in the field equations can eliminate the diver- 
gences that have cropped up in both the preceding problems. 

Despite its uncertainties, the analysis of Rosen's theory in this 
section leads one to the following tentative conclusion: if v^ < 1 , the 

speed of gravitational radiation is the ultimate speed for particles of 
nonzero rest mass — a "speed limit" enforced by the emission of gravita- 
tional radiation. Hence, observations of relativistic particles can be 
used to place limits on the speed of gravitational radiation (Eq. (l.3a)) 
and on the cosmological boundary values [Eq. (1.3b)]. 
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4. OBSERVATIONAL CONSTRAINTS 


Tlie hlghes t-energy particles In the vicinity of the ear^h are ultra- 

high-energy cosmic rays, which have been detected at energies exceeding 

20 

10 eV (see (22] for a review of the observations). At these very high 

energies, cosmic rays are not observed directly; rather, they are detected 

by the air shower they produce as they enter the atmosphere. The energy 

assigned to the primary particle In a given event is somewhat uncertain, 

since it is derived from a model for the shower. However, an energy of 
19 

3*10 eV seems reasonably firm. 

This energy estimate, even if correct, is not a measurement of veloc- 
ity. One obtains a velocity by using the familiar relation E ■ ym . Hov- 

o 

ever, one might expect this relation to fall in Rosen's theory, because a 

particle's gravitational binding energy might diverge as v approaches v . 

8 

Indeed, an analysis using the linear approximation suggests that the energy 
of a particle, as measured by an observer at rest in the LURF, diverges 
logari thmically : 

E - ym 4- >q log{2v g [l - (y/Y g ) 2 ] , (4.1) 

where is the gravitational binding energy when the particle is at rest 
in the LURF, and q is a dinu nsionlesu quantity which depends on the structure 
of the particle. This divergence is one more reason why particles cannot 
exceed the speed of gravitational radiation. 

The logarithmic divergence (4.1) is slow enough that it does not 
interfere with interpretation of the cosmic-ray observations, ff a particle's 
speed is so close to v^ that the binding-energy term in (4.1) dominates, 
then Eq. (3.6) predicts that the particle will radiate away almost all its 
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energy as gravitational radiation. It will not produce the observed shower 
of particles. 

A more serious uncertainty results from the Inability to identify the 

primary particle. The most likely candidates are protons or, perhaps, alpha 

particles; however, the possibility of heavier nuclei — perhaps nuclei near 

19 

iron — has not been ruled out. For a proton at 3*10 eV, the limit (1.3a) 
on the speed of gravitational radiation near the earth today is 


-22 

1 - v < 5 * iU 

g ' 

For an iron nucleus at the same energy, the limit is a bit weaker: 

• 1 H 

1 - v < 1 * 10 

K ' 


(4.2a) 


(4.2b) 


Since v^ Increases toward the galactic center, these limits also apply at 
any point closer to the galactic center than the earth. 

Equations (4.2) actually hold not only at the earth hut also in those 
regions traversed by the cosmic rays after their initial acceleration. 
Unfortunately, the point of origin of ultra-high-energy cosmic rays is un- 
certain. Their Larmor radii in the galactic magnetic field are much larger 
than the thickness of the galactic disk. This, together with the lack of 
anisotropy in the observed events [23], means that, *■ £ they are galactic in 
origin, they must come from a distance less than the thickness of the disk 
('V 200 pc). It seems more likely that they are ext ragalact ic , in which 
case Eqs. (4.2) probably apply out to a distance of at least 100 Mpc . 

Earth-based observations of relativistic particles also provide an 
upper bound on the value of (v ) today [Eq. (1.3b)]. This limit is 
considerably less stringent than the limit on v^ because it is determined 
by the Newtonian potential at the earth, which is dominated by the galactic 
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potential A 1 particle with y > 10 — a medium-energy coawic-rav 

proton or electron, a poaltron or electron produced In a high-energy colli- 
sion at Fermi Lab or CFRN or circulating In a storage ring at SLAC or DESY — 
yields the same limit: 


< 4U 


today 


gal 


^ 3 * 10 


-6 


(4.3) 


Here I have used a galactic mass of 1.4*10^M at a distance of 10 kpc. 

# 

For positive values of a,, this limit (valid only In Rosen's theory) is 
almost three orders of magnitude better than the best previous limit, ob- 
tained by searching for anomalous earth tides [24]. 

There la a possibility that the Newtonian potential of the Virgo clus- 
ter at the earth is as large as the galactic potential. However, there Is 
considerable uncertainty in estimating the mass of the Virgo cluster, and 
the two potentials are comparable only for the largest estimates. In anv 
case, Including the potential of the Virgo cluster is not llkelv to degn.de 
the limit (4.3) by more than e factor of two. 

Compact radio sources at substantial red shifts provide information 
about the speed of gravitational radiation in the past. They enlt a power- 
law radio spectrum which is thought to be Incoherent electron-synchrotron 
radiation; the spectrum has a low-frequency turnover attributed to synchro- 
tron self-absorption. The Lorentz factor of the electrons can be estimated 
from the brightness temperature T| at the turnover frequency: Y ^ (kT^/n.^), 

where in 1 the rest mass of an electron. Jones, O'Dell, and Stein [251 
e 

have developed a detailed model for compact, nonthermal sources, including 
the effects of synchrotron self-absorption am’, synchrotron self-Compton 
radiation. Burbidge, Jones, and O'Dell [26] have applied the model to 
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several compact sources, some of which have more than one component (see 
their Tables 1-3). For nine of the ten sources In their sample, they pro- 
vide (for one or more of the components) a red shift, an angular diameter 
determined by VLBI, a size determined from the angular diameter by placing 
the source at Its red-shift distance, and a Lorentz factor determined by 

the model. To estimate a Newtonian potential for each component, I have 

9 

assumed a mass of 10 M — a mass larger than or of order those usually thought 

o 

to be associated with active galactic nuclei; and 1 have assumed a constant 

(nonevolving) gravitational constant G. As an example, consider the source 

with the largest red shift in their sample — PKS 213^ + G0.U at z* 1.936. The 

-5 

estimates for one of its two components are y — 590 and U — 1 x 10 , which 

-5 

implies O o <Uxl0 [Eq. (1.3b)]. Similar considerations for the other 
sources provide upper bounds on v , «t a variety of red shifts; considering 
all these limits together, one can conclude that 

'r-, ~ 5 x 10* U for 0 ^ a 4 2 . (U.l») 

No other observation provides information about the value of (v^ c ) * n t ^ ie 
past. 

There are considerable uncertainties in estimating the Newtonian poten- 
tials which go into the limit (U.U). The masses and radii of the sources are 
uncertain; in addition, the gravitational "constant" does evolve in Rosen's 
theory, so that its value in the past depends on the cosmological model. 

Because of these uncertainties, the limit (U.b) has been chosen conservatively; 
with the above assumptions, all but one of the slot sources at z >0.1 provide 
a limit et least an order of magnitude stronger than ( 1* . U ) . 

5. CONCLUSION 

Rosen [27] has recently modified his "bimctric" theory. In the modi 
fied version, the "background" metric T1 is no longer required to be flat; 
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instead, It is required to bo a space of constant curvature. Cosmological 
models are affected by tit i m modification, but local gravitation physics 
Is not, except Insofar as It Is Influenced by cosmological boundary values. 
The analysis In Sections 2 and I remains the same, and the limits obtained 
In Section 4 apply to the new version of the theorv . 

The analysis In this paper has been restricted to Rosen's theory, but 
the results obtained are likely to have far rfider applicability. There are 
numerous metric theories of gravity which predict different speeds for 
gravitational radiation and light. Typically In such theories, the differ- 
ence in speed Is produced Just as In Rosen's theory: light propagates 

along "light" cones of the physical metric, while gravitational radiation 
propagates along "light" cones of a flat, "background" metric. In all such 
theories, one expects the speed of gravitational radiation to have a foim 

similar to that in Rosen's theory: v * v (l + 2£,U). where v is det« r- 

g gc «c 

mined by cosmological boundary values and f. Is a constant of order unit) . 
The important question Is whether emission of gravitational radiation re- 
stricts particles to speeds less than v^. Although detailed calculations 
are necessary in each theory, one can give a general argument, based on 
the analysis in Rosen's theory, for the existence of the "gravitational 
speed limit." 

Whenever a particle exceeds the speed of propagation of a "radiation" 
field to which it is coupled, one expects a shock wave to form. One can 
think of numerous examples, such as the shock wave produced by supersonic 
motion In an acoustic medium and electromagnetic Cherenkov radiation. In 
these familiar examples, the radiation does not diverge because the shoe* 
front Is not absolutely sharp; it is spread out over some length d charac- 
teristic of the medium through which the radiation Is propagating. Tills 




"blurring" of the shock front cuts off the radiation at frequencies > v v , avt / < ** 
In the gravitational case, the "medium" is spacetime itself, or more ac- 
curately, the "background" structure or. spacetime which determines the 
speed of gravitational radiation. The gravitational "medium" has no small- 
scale structure to blur the shock front. Thus, there is no high-f requen :y 
cutoff (unlLSs quantum mechanics Introduces one), and the radiation does 
diverge. 

This argument makes it seem quite likely that any theory with a 
variable speed of gravitational radiation must confront the limits obtained 
in Section 4. If so, Eqs. (4.3) and (4.4) can be used to constrain the 
cosmological boundary values in any such theory. (In general, the Qg of 
these limits is not a PPN parameter; it is simply a parameter related to v^ c 
by Eq. (1.2a),] In additi n, Eqs. (U.2) provide a general, theory-indep ndent 
lower bound on the speed of gravitational radiation near the earth. 
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FOOTNOTES 


Throughout I use the summation convention. Creek indices running from 0 
to 3 and Latin indices from 1 to 3. The signatures of the metrics are +2. 
A semicolon ( ;) denotes a covariant derivative with respect to g^ )t •> 
vertical bar (|) a covariant derivative with respect to n , and a comna 
an ordinary partial derivative. Units are chosen so that the speed of 
light c ■ 1 . 
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FIGURE CAPTION 




Figure 1. A "snapshot", taken at time t, of the Cherenkov cone produced by 
a particle moving with uniform velocity v ■ 2v^ along the z-axis. Thp 
particle Is at the apex of the cone. The angle 0 C between the normal 
to the cone and the z-axls is given by cos 0 » (v /v) . The shaded 
region is the pulse of Cherenkov radiation produced by a particle which 
radiates from t ■ 0 to t ■ T . 
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